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A new formula for the force vs extension relation is derived from the discrete version of the so 
called worm like chain model. This formula correctly fits some recent experimental data on polymer 
stretching and some numerical simulations with pairwise repulsive potentials. For a more realistic 
Lennard- Jones potential the agreement with simulations is found to be good when the temperature 
is above the 9 temperature. For lower temperatures a plateau emerges, as predicted by some recent 
experimental and theoretical results, and our formula gives good results only in the high force 
regime. We briefly discuss how other kinds of self-interactions are expected to affect the elasticity 
of the polymer. 



PACS numbers: 

I. INTRODUCTION 

The behaviour of a single polymeric molecule under 
stretching has become a popular subject among exper- 
imental physicists (see e.g. [1], for a review) and has 
attracted the attention of theoreticians who have intro- 
duced several models to explain these results. An exhaus- 
tive description on this topic is hard to achieve within a 
single model, since the experimental conditions strongly 
affect the results. So, the most common theoretical ap- 
proaches are usually only a first step approximation. 

Here, we focus our attention mainly on the pulling of a 
polymer under good solvent conditions [2] , in which case 
we give an analytical prediction for the force versus ex- 
tension curve. We also discuss other cases by means of 
numerical simulations and direct connections to experi- 
ments. Usually, the models adopted describe the poly- 
mer as an elastic (ideal) chain, where self-avoidance is 
not taken into account. Typically, freely jointed (FJC) 
or worm like chain (WLC) are studied [3]. 

The former describes a chain of beads connected by 
links of constant distance, whereas the latter introduces 
an intrinsic stiffness between two consecutive bonds and, 
in particular, is shown to correctly describe a wide range 
of experimental results on double-stranded (ds) DNA, 
single plasmid and lambda phage DNA [4-6]. In this 
case, the large force behaviour 1 — z/Lc ~ 1/VF is found, 
where z is the elongation along the direction of the force, 
Lc the contour length of the polymer and F is the applied 
force [4, 7]. Let us note that in these cases the contin- 
uous version of the WLC model is always used, where 
the persistence length [4] Lp is very large, compared to 
the base separation (roughly one persistence length is 150 
base pairs). 

Here, we focus on the discrete version of the WLC 
model, that satisfactorily describes the pulling behaviour 



of a polymer in good solvent. We broadly identify three 
regimes in the force vs extension curves obtained in our 
analytical and numerical calculations. The low force (or 
low stretch) regime is highly affected by the details of 
the interactions between the beads (which in nature are, 
e.g., caused by the different concentration of ions in so- 
lution). This regime is discussed only marginally here, 
as, though it is potentially very interesting, experimen- 
tal data in this range of forces are quite rare and not 
precise enough to allow a comparison with theoretical 
predictions. There is then a second regime, of interme- 
diate stretches or forces, in which the force versus exten- 
sion characteristic curves obey approximately the laws 
predicted a few years ago in Ref. [4] by means of a con- 
tinuum theory of the WLC. Finally, for very large forces 
(beyond a polymer dependent crossover value), we get 
a universal, model-independent, freely-jointed-chain like 
behaviour. In Section II, we introduce the discrete WLC 
model and show how it can be related to the well known 
continuous version. In Section HI, we show a compar- 
ison between our theory and some recent experiments. 
Then, in Sections IV and V, we introduce some extension 
vs force curves obtained from Monte-Carlo and Molecu- 
lar Dynamics calculations for models of stiff and flexible 
polymers with a chosen potential between non nearest 
neighbour beads. We shall consider two different cases: 
a pure repulsive potential and a more realistic Lennard- 
Jones potential and we show that in the latter case our 
formula agrees with the numerical data only in the case of 
high temperatures or, for lower temperatures in the high 
force regime. Below the point a force plateau appears, 
as predicted by some recent theoretical works [8-11] and 
confirmed by experiments [5]. Finally, Sections VI and 
VII are left for discussion and conclusions, respectively. 



2 



II. THE MODEL 

Our model describes a chain of beads, where the dis- 
tance between the nearest neighbours is kept fixed (we 
can put it equal to 6) and a suitable stiffness K is intro- 
duced. Then, the corresponding Boltzmann weight reads 
[12-14]: 
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where /3 = 1/T, T being the temperature in units of 
Boltzmann constant, tj = Fj — rj_i (rj being the position 
vector for the i-th bead, i = 1,...,N and A'' is the total 
number of beads) and F is the applied force which defines 
the z-direction. The partition function for the model 
described by (1) is 
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The average elongation {z)n along the stretching direc- 
tion is 
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In the thermodynamic limit (A^ 
behaviour of Eq. (3) gives 



C = lim 



N 



Nb 



= 1 - 



(3) 

oo), the large force 
(4) 



The continuum approximation of Eq. (4) is obtained 
with the following substitutions [15]: 



PK ^ Lp/b, 



(5) 



in the formal limit b 0, where Lp is the persis- 
tence length. The final result ( = 1 — \/2^J (iLpF does 
agree with the celebrated result of Marko and Siggia [4] . 
However, our result is more general, since it predicts a 
crossover force Fc = 



4Lp 

WLC-Uke behaviour: 1 - C ~ l|^/F, F <^ 
FJC-like behaviour: l-C^l/F, F ^ Fc 



(6) 



Let us notice that the validity of the continuum approx- 
imation proposed in Rcf. [4] is not simply related to the 
value of the dimensionless ratio b/Lp but rather to F/F^- 
From Eq. (4), we deduce that: 
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The low force behaviour of Eq. (3) is 

l3bF 1 + v{Lp/b) 



3 l-y{Lp/b)' 



(7) 
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where y{x) = coth(x) — l/x [16]. Again, in the limit 
6^0, Eq. (8) agrees with the result of Marko and 
Siggia [4]. Matching Eqs. (7) and (8), we obtain the 
following interpolation formula: 
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It is easy to verify that Eq. (9) correctly reproduces the 
right large and small force behaviours, Eqs. (7) and (8) 
and that in the continuum limit 6 ^ we obtain the well 
known interpolation formula [4] : 
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In next Section, we shall apply our formula, Eq. (9), to 
two recent experiments. The first discusses the stretch- 
ing of a single plasmid DNA molecules, to which the for- 
mula (10) was previously applied with success [5]. Nev- 
ertheless, our formula predicts a non trivial value for the 
parameter b, that gives an estimate for the intra bead 
distance. Then, the second experiment [17, 18] demon- 
strates that Eq. (9) gives the right large force behaviour. 

In the following, we shall define C, = zjL^, where z is 
the elongation along the direction of the force and = 
Nb is the contour length of the polymer (see also Eq. 
(4)). 



III. COMPARISON WITH TWO RECENT 
EXPERIMENTS 

Let us consider the experimental data reported in the 
plot at the top of Fig. 3 of Ref. [5]. The authors 
considered the elastic response of a single plasmid of 
DNA molecules, probed using optical tweezers [19]. They 
found that, according to the environmental conditions, 
it can be very different. In particular, for condensed 
molecules, the stretching curves display a stick-release 
pattern, where the DNA molecule can be described as 
a succession of different WLC's of different contour and 
persistence lengths. We are mainly interested to the case 
of uncondensed molecules, whose stretching pattern is 
reproduced in Fig. 1. 

As described in the caption, the contour and persis- 
tence lengths obtained with our formula, Eq. (9), and 
with Eq. (10) are perfectly compatible. Nevertheless Eq. 
(9) predicts a non trivial value for the intra bead effec- 
tive distance 6 = 2.5 nm, that corresponds to 7 — 8 base 
pairs. Let us notice that this matches the DNA hydration 
thickness (here, we have used ^ 0.34 nm as the distance 
between two consecutive base pairs [5]). 

If we take T = 300 K, the crossover force F^ ~ 100 
pN and F/F^ < 0.16 for the data of Ref. [5], thus jus- 
tifying the use of the interpolation formula, Eq. (10). 
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FIG. 1: (+): experimental data taken from Ref. [5]. Long 
dashed line: fit with the curve, Eq. (9), yielding Lc = 1324 

nm, Lp = 38 nm and b — 2.5 nm. Short dashed lino: fit with 
the curve, Eq. (10), yielding Lc = 1324 nm and Lp = 38 nm, 
see Ref. [5]. 




FIG. 2: (+): experimental data taken from Ref. [18]. Con- 
tinuous line: fit with the curve, Eq. (9), yielding Lc = 2.31 
(in units of the contour length l'^^ of the equivalent dsDNA 
molecule observed in 10 mM PB [18]), Lp = 0.26 nm and 
6 = 0.12 nm. 



However, as pointed out in Eq. (6), the discrete nature 
of the chain does emerge, when F > F^. Notice also that 
for forces considerably smaller than F^ double-stranded 
DNA would undergo an overstretching transition [20], 
where a more sophisticated theory is needed [3] . 

An interesting question is how this treatment may be 
applied to single-stranded DNA (ssDNA). On one hand, 
if we keep as physical parameters the persistence length 
of ssDNA ~ 1 nm. and as the equivalent of h the sepa- 
ration between two phosphates, i.e. 0.5 nm roughly, we 
would end up with a crossover force again of the order of 
70 pN. Data in this regime do exist [17, 18], and suggest 
that the WLC grossly fails to fit the data [18, 21]. In 
fact (see Fig. 4, Ref. [18]), the authors pointed out that 
the corresponding fit with Eq. 10 gives good results in 
the large force regime, but with a calculated persistence 
length (~ 0.21 nm) which is clearly not physical. Our 
equation does not do much better for low and intermedi- 
ate force, in which case, as shown in Ref. [18], evidently 
the self-interactions dominate the behaviour. Still a large 
force fit, even if a bit dependent on the contour length 
which we choose, suggests that the large force exponent 
in the (log(l — C), log(F)) plane is —1 as predicted by our 
model (see Fig. 2). The calculated fitting parameters Lp 
and h (see the caption of Fig. 2) give F^ 300 pN, which 
is consistent with our Eq. (9) (see Fig. 2). 

In next Section we introduce some Monte-Carlo calcu- 
lations and compare them to Eqs. (9) and (10). 



IV. MONTE-CARLO CALCULATIONS 
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FIG. 3: (+): Monte-Carlo data for a stiff chain oi N = 100 
beads and stiffness K = 40, for /3 = 1. Long dashed line: Eq. 
(9) with Lp = 40 and fe = 1. Short dashed Une: Eq. (10) with 
Lp = 40. The crossover force Fc, see Eq. (6), is also shown. 



and the curves given by Eqs. (9) and (10) (long and 
short dashed line, respectively), for Lp = 40, and f3 = 
1. As observed, the agreement is perfect only for Eq. 
(9). In fact, the discrete nature of the chain emerges 
around F = F,. ^ 160 (in the chosen units) and the 
WLC approximation is no more valid. 

To render the model more realistic, let us introduce a 
short range repulsive potential Vrep{r) between non con- 
secutive nearest neighbour beads. We have adopted the 
following functional form: 



As already said, our model is a stijf chain described 
by the Boltzmann weight, Eq. (1), where the intra bead 
distance b is now kept fixed to 1. 

In Fig. 3 we have plotted the Monte-Carlo data (-I-) 



Vrep{r) = 



.12 



(11) 



where r is the distance between two beads. Then, we 
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FIG. 4: (+): Monte-Carlo data for a stiff chain oi N = 100 
beads, stiffness K = 30 and interaction potential given by Eq. 
(11), for /3 — 1. Long dashed line: fit obtained from Eq. (9) 
with Lp ~ 31 (6 = 1). Short dashed line: fit obtained from 
Eq. (10) with Lp ~ 34.9. 



have to multiply the Boltzmann weight, Eq. (1), by 

where we have defined = | r j — r ^ | . 

In Fig. 4 we have plotted the Monte-Carlo results (+) 
for K = 30, together with the two fitting lines obtained 
from Eqs. (9) and (10) (long and short dashed lines, 
respectively). The corresponding persistence lengths are 
Lp ~ 31 {b is kept fixed to 1) and Lp ~ 34.9. Again, 
we can observe that our formula works better than Eq. 
(10). Moreover, this last result is in agreement with some 
theoretical works [22, 23], that predicts that the net effect 
of the repulsive potential is to renormalize the persistence 
length, making it larger than the bare one. 

Now, let us introduce a more realistic potential, adding 
to the repulsive core an attractive part too. It is inter- 
esting to notice that the approaches of Refs. [22, 23] can 
not be generalized, since the attractive part introduces 
some instabilities and the perturbation scheme discussed 
there is no more valid. 

We have chosen the following Lennard- Jones kind func- 
tional form VLj{r): 

V,j{r) = Vo - ^) . (13) 

The parameters Vb and a are chosen in such a way that 
the minimum of the energy is located at r = rmin = 1-5 
and VLjir = Tmin) = —1- We expect that Eq. (10) does 
not work well in this case. Our goal is to test our formula, 
Eq. (9). 

Firstly, we show the effects of increasing stiffness on 
the force vs extension curves. To fix the ideas, we begin 
to consider a sufficiently high temperature, i.e. above the 
point [2, 16], whose location, however, is not exactly 



FIG. 5: (+): Monte-Carlo data for a stiff ch.am of iV = 100 
beads, stiffness K = IQ and interaction potential given by Eq. 
(13), for 13 = 0.3. Long dashed line: fit obtained from Eq. (9) 
with Lp ~ 2.12, b ~ 1.26. Short dashed line: fit obtained 
from Eq. (10) with Lp ~ 2.41. 

known. In Rcf. [24], the authors were able to numeri- 
cally determine the 9 temperature for a model of ho- 
mopolymer with a square well potential, whose depth is 
fixed to —1. They found that Tg ~ 3. Our Te should be 
somewhat smaller due to the stiffness. 

Initially, we fix /3 = 0.3 (T = 3.333...), with K = 10. 
The numerical results and the corresponding fitting line 
are plotted in Fig. 5. The corresponding fitting parame- 
ters are described in the caption. In this case we allowed b 
to be a free parameter and we have determined it through 
a best fit of the data even though in the simulated model 
6=1. The fit with our formula is surprisingly good, in 
contrast with the result obtained with Eq. (10). More- 
over, the predicted value for b is of the correct order. 
This means that the effect of the stiffness compensates 
the attraction due to the potential and the behaviour is 
similar to a FJC. Let us stress on the fact that if we fix 
fe = 1, the corresponding fits are considerably less precise. 

We have also simulated the case with K = 80. In Fig. 
6 we have plotted the numerical data together with the 
two fitting lines (the fitting parameters are reported in 
the caption). The agreement with our formula is again 
perfect, in contrast with Eq. (10). 

For both situations plotted in Figs. 5 and 6, we note 
that the main effect due to the potential is a considerable 
reduction in the persistence length due to the attractive 
part, which evidently renormalize also the intra bead dis- 
tance. 

As a check, let us notice that the correlation function 
(t, • t,+r) for the WLC is [16]: 

(t, • u+r) - vmr, (14) 

where the function y{x) has been defined above. In the 
continuum limit 

(ti • ti+r) = exp(-r/Lp). (15) 
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FIG. 6: (+): Monte-Carlo data for a stiff chain oi N = 100 
beads, stiffness K = 80 and interaction potential given by Eq. 
(13), for P = 0.3. Long dashed lino: fit obtained from Eq. (9) 
with Lp ~ 23.33 and b = 0.86. Short dashed line: fit obtained 
from Eq. (10) with Lp ~ 25.95. 
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FIG. 7: Plot of the correlation functions (ti ■ ti+r) vs r, cal- 
culated from the data of Figs. 3 to 6 and the corresponding 
theoretical predictions (continuous lines), given by Eq. (15). 



This result is exact for all r but holds in a more general 
context, like the case with interaction, in the large r limit. 
In Fig. 7 we have plotted the correlation functions (ti • 
ti+r) vs r for the data of Figs. 3 to 6 and compared 
them to the theoretical ones as given by Eq. (15). The 
corresponding persistence lengths are those found with 
our formula. The agreement is perfect. 

Let us now consider the case of /? = 0.5 (T = 2). In 
Fig. 8 we have plotted the case for K = 10. As it can 
be seen a first order phase transition emerges [9, 25]. In 
this case our formula is able to describe only the part of 
the plot at relatively large forces. 

A complete characterization of the force vs extension 
behaviour in this case needs a more complete theory 
and is beyond the scopes of this work. We reserve to 
study this important issue in a future work (see also Ref. 



FIG. 8: {+): Monte-Carlo data for a siiif chain of iV = 100 
beads, stiffness K = 10 and interaction potential given by Eq. 

(13), for l3 = 0.5. Long dashed line: fit obtained from Eq. (9) 
with Lp ~ 4.42 and b = 1.03. Short dashed line: fit obtained 
from Eq. (10) with Lp ~ 4.85. Let us remark that we have 
tried to fit only the range of data at large forces. 

[26] for some recent Monte-Carlo results about stiff poly- 
mers). 

V. MOLECULAR DYNAMICS CALCULATIONS 

In the Monte-Carlo simulations presented in previous 
Section we were working in fixed force ensembles, i.e. the 
force was fixed in each simulation and the averaged ex- 
tension of the chain was computed. It is interesting to 
check if in fixed stretch ensembles, i.e. when one fixes 
the extension and measures the averaged force, the re- 
sults remain unaltered. In the following we will present 
and discuss Molecular Dynamics (MD) simulation results 
in which we can perform a quasi-fixed stretch treatment. 
In these simulations (similar to what has been done for 
protein- like polymers [27, 28]) the two ends of the poly- 
mer are tethered by springs to two fixed points in space 
so its end-to-end vector can vary only little. The two 
fixed points are chosen such that a line connecting them 
is along the ^;-direction. We will consider two cases: a) 
when there are repulsions between two non-consecutive 
beads and b) when the beads interact via a Lennard- 
Jones (LJ) potential. Note that we do not include the 
intrinsic stiffness into the models considered in this Sec- 
tion (i.e. K = 0). 

Our homopolymers are modeled as chains of beads con- 
nected by springs. The harmonic potential for the chain 
connectivity reads 

Vi,i+i = ^k{n,i+, - bf , (16) 

where rj^j+i = [r, — rj_|_i| is the distance between two 

consecutive beads, b is fixed to 1 and k = 1444e is the 
chosen spring constant (e is a unit of energy). We use 
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the same spring constant k for the potential that tethers 
the ending beads to the fixed points. For the case of 
repulsion the following potential is used 



(17) 



where a = 1.32. For the case of attraction the LJ poten- 
tial takes the form 



/ \12 













(18) 



The model is studied by a MD method in which coupling 
to the heat bath is provided by the Langcvin equation 
(see e.g. Ref. [27] for details). In the following tem- 
perature will be given in units of e/ks where A;b is the 
Boltzmann constant. We consider chains of 60 residues 
and simulations are done for various temperatures. The 
averaged force (F) exerted on each of the fixed points and 
the averaged z-component of the end-to-end vector of the 
chain (z) are measured under equilibrium condition at a 
constant temperature. As we vary the position of one 
fixed point versus the other we are able to obtain {F) vs 
(z) under a quasi-fixed stretch condition (the fluctuation 
in z is considerably small). 

It should be noted that due to the springs used for the 
chain connectivity the contour length Lc is not a constant 
as one varies the stretching distance. We find that as the 
polymer approaches its full extension Lc may increase up 
to 5-10% depending on the temperature. We find that 
however if the relative extension ( is defined as {z)/{Lc), 
where (Lc) is the averaged contour length computed from 
the simulations at a given stretch z, then the {F) vs ( 
acquires very little changes as we increase fc by a factor 
of 2 or 3. 

In Fig. 9 we show the stretching data for the chain 
with repulsive potential given by Eq. (17) at T = 1 (at 
some other temperatures we found that the results are 
qualitatively similar). One can see that the data can be 
fitted quite well by Eq. (9) at high extension, but not 
at low extension. The least-squared fit to the data gives 
b ~ 1.05 and Lp ~ 1.04. Note that the value of b ob- 
tained from the fit is reasonably good though it is some- 
what larger than the bond length when no tension is ap- 
plied. This increment (of about 5%) however is in a good 
agreement with the increase in Lc under stretching and 
can be understood as due to the softness of the bond's 
potential. Thus Eq. (9) still gives the right physics. 
One can also notice that our homopolymer with repul- 
sive potentials corresponds to a worm-like chain with a 
low stiffness (Lp ~ b). We have computed directly the 
persistence length from the correlation (t, • U+r) by per- 
forming simulations for a free chain and found that Lp 
is approximately 2.35 at T = 1. The difference (of a fac- 
tor of 2) between the measured persistence length and 
the one obtained from the fit is due to the fact that the 
fit was good only at a high extension. In this regime 
the excluded volume effects are much smaller than in a 
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FIG. 9: Force (F) vs relative extension ^ for a 60-bead chain 
with repulsion studied by Molecular Dynamics simulations 
under quasi-fixed stretch conditions (see text) . is defined as 
{z)/{Lc) where (z) and (Lc) are the averaged extension and 
the averaged contour length obtained in the simulation. The 
data points shown are for temperature T = 1. The fit to Eq. 
(9) (solid line) yields b ~ 1.05 and Lp ~ 1.04. Eq. (10) fails 
to fit the data (dashed line). 
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FIG. 10: Force (F) vs relative extension ^ for chain with 

Lennard- Jones (LJ) potentials under quasi-fixed stretch con- 
ditions. The fits of the data to Eq. (9) (solid lines) yield 
b ~ 1.01 and Lp ~ 0.88 for T = 4, and 6 ~ 1.0 and Lp ~ 0.5 
for T = 2. Eq. (10) fails to fit the data at both temperatures 
(dashed lines). 



free chain, thus resulting in a lower persistence length. 
It should be noted that Marko and Siggia's formula (Eq. 
(10)) completely fails to fit the data. This is related to 
the fact that for chains with low stiffness the continuum 
limit 6 — » cannot be applied and the discreteness of the 
chain should be taken into account. 

We turn now to discuss our results for homopolymcrs 
with attraction. Fig. 10 shows the force vs relative ex- 
tension for a 60-bead chain with the Lennard-Jones po- 
tential given by Eq. (18). We find that Eq. (9) can 
reasonably fit the data for temperatures higher than the 
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collapse transition temperature Tg whereas at lower tem- 
peratures it can fit the data only at high extension. The 
Tq temperature is estimated to be between 2 and 3 for 
this system. As shown in Fig. 10 at T = 4 Eq. (9) 
appears to fit both high and low extension data yield- 
ing b ~ 1.01 and Lp ~ 0.88. At T = 2 the fit can be 
given only at large extension, which yields 6 ~ 1.0 and 
~ 0.5. Note that at both temperatures the value of b 
obtained from the fits is quite consistent with the correct 
bond length. Like for the case of repulsion, due to the 
low stiffness induced by the potentials, the Marko and 
Siggia formula does not fit the data as our new formula 
does. 

It should be noted that when the temperature is about 

the collapse transition temperature Tg we observe a 
plateau in the force vs extension (see the case of T = 2 
in Fig. 10), similar to what seen in the previous Section. 
However, at temperatures lower than Tg, in the fixed 
stretch ensemble the force vs extension curve shows the 
typical behaviour with hysteresis phenomenon [27, 28]. 



VI. DISCUSSION 

Though our main results refer to a situation in which 
attractive interactions are not very important, it is useful 
to further comment on how our results change in the 
presence of such a situation in a real experiment. Such 
attractions change the picture depicted so far as, if they 
are strong enough, they can cause a phase transition in 
the molecule. 

An effective self-attraction, like the one considered here 
by means of a 6 — 12 Lennard- Jones potential, arises due 
to hydrophobic interactions in polypeptides and due to 
a suitably large concentration of polyvalent counterions 
in the case of double-stranded DNA. If this is the case, 
the polymer is in poor solvent conditions or equivalently 
below the point (in the nomenclature of the models 
discussed above). In these cases at zero stretch the force 
attains (in a thermodynamically long chain) a non-zero 
value followed by a force plateau for long molecules [9]. 
This indicates the presence of a first order transition. 
We have already seen in the previous section that in this 
situation (i.e. in poor solvent conditions) the two ensem- 
bles, fixed force and fixed stretch, are not equivalent and 
indeed in the experiments the stretching curves of small 
DNA's or of DNA's in presence of a high concentration of 
polyvalent counterions in reality show peaks. A thorough 
explanation of this effect includes polyelectrolyte model- 
ing, finite size corrections as well as dynamical effects and 
can be found in Refs. [9, 10]. 

A different scenario holds for ssDNA and for RNA (this 
scenario would also be retraced if a self-attractive poly- 
mer such as dsDNA in presence of condensing agents 
is restricted conformationally to stay in a quasi-two- 
dimensional film). In these cases the single molecule 
phase transition is second order and thus at low stretch 
the characteristic force curves are non-zero and then rise 



smoothly. Furthermore, the attraction in ssDNA and 
RNA is brought forth by the base pairing interactions 
between bases far apart in the chain. If sequence disor- 
der is neglected, the low force regime can be written as 
follows: 

C^(,/-.fc)^/^-^ (19) 

where fc is the critical force. The exponent A is between 
0.5 (the value of native branched polymer-like configura- 
tions arising from base stacking [29]) and 1 (for a polymer 
in a good solvent). While present day experiments show 
that as the extension goes to the force is non-zero [30], 
the data are not clean enough to allow discriminating be- 
tween the exponents above. To be noted that within a 
simplified theory [30], it was possible to find a law im- 
plying A = 1/2. 

VII. CONCLUSIONS 

Here, wc have revised the well known WLC model that 
correctly describes the behaviour under pulling of a stiff 
polymer. We have pointed out that its discrete version 
has a different large force behaviour respect to the contin- 
uous version by Marko and Siggia [4] , Eq. (6). We predict 
a crossover force between these two different regimes. It 
should be noticed that a recent paper by Livadaru et 
al. [31] reports a similar result. However, here we have 
given a simpler formula, Eq. (9), which can be tested on 
real polymers as well as on numerical simulations, where 
self-interactions are present. 

So, firstly, we have used it to fit some experimental 
data on dsDNA, Fig. 1, and ssDNA, Fig. 2. In the 
former case, where Eq. (10) was already applied with 
success, we find that Eq. (9) predicts a not trivial in- 
tra bead distance and a crossover force, which is much 
greater than the experimental forces. This justifies the 
approach by Marko and Siggia. In the latter case, we 
observe a clear crossover between the two regimes. In 
fact, now the crossover force is small and the assumption 
leading to Eq. (10) is no longer justified. 

Then, we have performed some Monte-Carlo simula- 
tions to verify the validity of our formula. We have an- 
alyzed different kinds of intra bead potentials. For a 
short range repulsive potential we have correctly found 
that our formula is in good agreement with the numeri- 
cal data, predicting a renormalized stiffness [22, 23]. For 
a more realistic Lennard-Jones potential the situation is 
more complicated. For temperatures above the 9 point, 
Eq. (9) gives a good fit. The main result is that, now, 
both the persistence length and the intra bead distance 
b are renormalized by the potential. For temperatures 
below the 9 point, our formula agrees with the numerical 
data only at large forces. 

We have performed also simulations to examine the va- 
lidity of our new formula for the pulling of flexible poly- 
mers under quasi-fixed stretch conditions. This kind of 
simulations has been done by using Molecular Dynamics 
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methods. Like in Montc-Carlo simulations we considered 
chains with repulsive and attractive Lennard- Jones po- 
tentials but now the chains have no intrinsic stiffness. It 
is found that for both types of interaction our new for- 
mula works very well at high extension and it correctly 
predicts the intra bead distance. The case of attraction 
is even more interesting since when the temperature is 
above the 6 point both the the high and low extension 
stretching data can be fitted reasonably well to Eq. (9). 
It is shown that the Marko and Siggia formula does not 
fit the numerical data for the chains considered. The 
reason for this is that the modest stiffness induced by 
the Lennard- Jones potentials is not sufficiently large to 
yield the continuous WLC behaviour as described by Eq. 



(10). The results indicate that Eq. (9) can be used to 
characterize elastic behaviours of a much wider range of 
biomolecules. 

Finally, we have discussed how an attractive potential 
can modify the pulling behaviour of a polymer under poor 
solvent conditions. 
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